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The koala Phascolarctos cinereus is the only member of the once diverse marsupial family Phascolarctidae 
to have survived the Last Glacial Maximum. A climate envelope model for P. cinereus was developed 
to predict the range for this species at present and at the Last Glacial Maximum. The model was 
compared to the contemporary koala records and the known fossil records of P. cinereus during the 
Quaternary.The predicted current core range for koalas was concentrated in southeast Queensland, 
eastern New South Wales and eastern Victoria. At the Last Glacial Maximum their predicted core 
range contracted significantly to southeast Queensland and northeast New South Wales. Our findings 
concord with other studies that find species experienced range contractions during glacial maxima. 
In the context of the future conservation planning for koalas in the wild, our historical perspective 
demonstrates the past adaptations of koalas to changes in climate and their probable range contraction 
to climatic refugia.The future survival of wide-ranging specialist species, such as the koala, may depend 
on identifying and protecting, future climatic refugia. 
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Introduction 

There has been a longstanding debate regarding the 
use of the theory of glacial refugia to describe current 
patterns of biodiversity (Willis and Whittaker 2000), 
but recently there is an increased interest in the 
importance of refugia in maintaining diversity under 
changing climates (Sommer and Zachos 2009). From 
an evolutionary perspective, the term ‘refugia’ refers to 
core habitats where organisms are able to persist during 
periods in which their wider geographic distribution 
becomes uninhabitable (Svenning et al. 2008; Fleikkinen 
et al. 2009). Flowever, this is a biologically complex and 
debated concept (e.g. Klicka and Zink 1997; Knapp and 
Mallet 2003). Most living organisms have experienced 
significant changes in their distributions due to global 
climatic fluctuations and particularly due to major ice 
ages over the past three millions years (Flewitt 2000). 
The Last Glacial Maximum (LGM) was estimated to 
have been at its maximum from at least 22-19 thousand 
years BP (Yokoyama et al. 2000). In Australia, ice-sheet 
cover was not extensive during the LGM and there is 
evidence, for example, that Eucalyptus regnans, was able 
to persist either locally during the harsh environmental 
conditions or by contracting to near-coastal refugia 
(Nevill et al. 2009). The identification of refugia for plant 
and animal species may have important implications for 
the conservation of biodiversity under climate change. 


The modern koala Phascolarctos cinereus is an arboreal 
folivorous marsupial endemic to Australia. It currently 
occupies a widely distributed but fragmented range 
in eastern Australia (Fig. 1). The evolution of the 
koala family (Phascolarctidae) was driven by an 
adaptation to open sclerophyll forests as the continent 
became drier (Hocknull et al. 2007) and to feeding 
almost exclusively on the leaves of eucalypt trees. This 
increased dependence on eucalypts through time is 
demonstrated for example, in changes to koala skulls 
(Louys et al. 2009). The fossil records of Phascolarctidae 
reveal a unique dentition of double W-shaped blades 
on each molar for slicing leaf matter. There are subtle 
variations in the koala’s dental nomenclature which 
enables taxonomic differentiation (Long et al. 2002), 
but the dental morphology of the modern koala is 
only marginally more specialised than its ancestors 
(Archer et al. 2000). The fossil records show that the 
distribution of koalas at times in the distant past 
was far wider than at present (Fig. 2) and also that P. 
cinereus is the last survivor of this once diverse family 
containing six genera and at least 18 species (e.g. 
Stirton 1957; Merrilees 1967; Archer 1972; Marshall 
1973; Bartholomai 1977; Pledge 1992; Black 1999; 
Piper 2007; Price et al. 2009). 
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The role of humans as well as climate in determining the 
biogeography, evolution and extinction of faunas is an 
area of ongoing research, in particular the investigation 
of fauna-rich fossil caves to ascertain faunal responses 
to climate change during the Pleistocene (e.g. Hocknull 
2007; Prideaux et al. 2007). These authors have revealed 
significant population fluctuations driven by glacial- 
interglacial cycling, but stable species composition 
for 500 thousand years, with the exception of the 
megafauna that went extinct around 70 - 50 thousand 
years ago (Long et al. 2002). Research suggests that 
climate change alone was unlikely to have been the 
principal cause of the extinction of the megafauna and 
although Aboriginal hunting may have been a factor 
in determining the distribution and abundance of P. 
cinereus, there is no doubt that populations have been 
significantly affected by the European occupiers of 
Australia for the fur trade in more recent times (e.g. 
Fowler 1993; Hrdina and Gordon 2004). 

The modern koala’s primary food source is a regionally - 
variable selection of, primarily, Eucalyptus and Corymbia 
tree species (e.g. Moore and Foley 2000; Ellis et al. 
2002). Koala populations are currently under threat 
from the synergistic threats of habitat loss, disease, 
car collisions and dog predation in coastal regions of 
Queensland and New South Wales (Jackson et al. 1999; 
Dique et al. 2003; McAlpine et al. 2006; Rhodes et al. 
2008; Lunney et al. 2009). For example, the abundance 
of koalas on the Koala Coast of southeast Queensland is 
estimated to have declined by 64% in a ten year period 
(1999-2008) (State of Queensland 2009). In western 
semi-arid and arid regions, koalas have been shown to 
experience significant local population declines during 
drought and heatwaves (Gordon et al. 1988). 

Climatic refugia may have played a key role in 
the survival of the extant koala species during the 
extreme climates of the LGM. Here we examine this 
hypothesis through a review of the past occurrence of 
P. cinereus, and comparison of its present and predicted 
LGM distribution using bioclimatic modelling. We 
discuss the implications for the management and 
conservation of contemporary koala populations 
under future climate change. 

Methods 

Approach and assumptions 

Bioclimatic models are useful tools for predicting the 
responses of organisms to climatic changes because 
they focus on the climatic thresholds that limit species’ 
distributions in past and present time (Mesquita and 
Sousa 2009). However, it must be remembered that 
models are simply estimates of species’ potential 
ranges. We used Bioclim from the ANUCLIM software 
package (Houlder et al. 1999) to develop a species 
profile or ‘climate envelope’ of the distribution of 
the modern koala. Bioclim has been used in a large 
number of studies world-wide (e.g. Manning et al. 2005; 
McKenney et al. 2007; Green et al. 2008), including 
palynological studies (e.g. Moss and Kershaw 2000; 
Gallagher et al. 2003). As a climate envelope model. 


Bioclim is independent of other variables that can 
influence koala presence such as habitat loss and 
urbanisation. Bioclim has been found, when compared 
to other climate envelope modelling methods, to be 
less prone to overprediction and therefore useful 
as a conservative approach (Hijmans and Graham 
2006). In another comparative modelling study on 
18 terrestrial taxa (insects, amphibians, reptiles, birds 
and mammals), Bioclim attained increasingly high 
concordance with other modelling methods, with 
increasing sample sizes (Hernandez et al. 2006). For 
this study we used a large sample size (n = 3,332). 

Due to the broad temporal span of the P. cinereus fossil 
records (< 1 to 569 thousand years ago) (Price 2008), 
these fossil data were not considered reliable for use in 
our modelling. We therefore used the contemporary 
distribution data of P. cinereus with its contemporary 
climate envelope and then used it to predict the past 
distribution of this species at the LGM. Determining 
regional climates at the LGM is a complex task. 
For example, Williams et al. (2009) found that mean 
annual temperatures in tropical northern Australia 
ranged from 3 to 7 ° C below present values and up 
to 8 ° C lower in mainland southeast Australia, with 
winter rainfall reduced by at least 30 % in tropical 
northern Australia. We therefore assumed generalised 
continental-scale decreases in rainfall and temperature 
for the study area, which may not reflect regional 
variations that occurred during the LGM. This study 
therefore examines potential macrorefugia (Ashcroft 
2010) at the LGM. 

Study area and populations 

We used contemporary koala distribution data from 
databases held by the States of Queensland, New 
South Wales and Victoria (Fig. 1). Koalas currently also 
inhabit parts of southeast South Australia (Fig. 1) where 
they have been re-introduced after virtual extinction 
(Jackson 2007), an enclosed reserve in Western 
Australia, and various islands around the south and 
east of the continent. To avoid uncertainty regarding 
whether such populations have been translocated or 
are native, we excluded these populations from the 
analysis and only considered areas containing natural 
populations. We defined ‘natural’ as populations 
that are not derived from solely translocated and/or 
introduced individuals. 

Data and climate envelope development 

Bioclim uses mathematically derived climate surfaces, 
estimated from mean monthly climate estimates of a 
network of meteorological stations across Australia, to 
develop the climate envelopes based on the latitude, 
longitude and elevation. We used the modern koala 
sighting records over the past 100 years. The koala 
sightings records were obtained from the Queensland 
Department of Environment and Resource 
Management, the New South Wales Department of 
Environment and Climate Change and the Victorian 
Department of Sustainability and Environment. 
Duplicate records, as well as obvious outliers indicating 
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Figure I. Modern koala range shown in black Study area shown in grey. Koala records sourced from: Queensland 
Department of Environment and Resource Management, New South Wales Department of Environment and Climate 
Change,Victorian Department of Sustainability and Environment and South Australian Department of Environment and 
Natural Resources. 



Figure 2. Published fossil sites containing relatively well-dated records of koalas since the Late Miocene, 11-5.3 million 
years ago (following Price 2008). 
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data entry mistakes, were removed manually by 
viewing the koala data points spatially. To reduce the 
possibility of geographical bias, such as a high density 
of records collected in highly populated areas, records 
were filtered by overlaying the occurrence records with 
a 10 kilometre grid and then one record was randomly 
selected from each grid cell (McKenney et al. 2007). 
This reduced the total number of koala records for the 
bioclimatic modelling from 32,029 to 3,332. 

Predictions: current range of P. 
cinereus 

We used Biomap, a component of Bioclim, to model the 
predicted current koala geographic range based on the 
climate envelope generated by Bioclim (temperature 
range -4 ° C to 37.7 ° C and annual precipitation 234 
to 2480 nun) (Table 1). For these predictions, we 
calculated the entire range (1-100 th percentiles) and 
the core range (10 th - 90 th percentiles). The entire range 
indicates locations that have bioclimatic parameters 
that fall within the maximum and minimum values 
of the climate envelope. This provides an estimate of 
potential koala habitats, including outliers and some 
may in reality not provide suitable habitat. The core 
range indicates the most climatically suitable locations, 
with points that have bioclimatic parameters that fall 
between the 10 th and the 90 th percentiles of the koala’s 
climate envelope. Flowever, the core range may include 
areas that, although optimal have, for example, been 
cleared for urban development or agriculture. 

Although Bioclim can utilise up to 35 climatic 
parameters, is has been suggested that using all 
of these can lead to over-fitting (Beaumont et al. 
2005). Following recommendations made by these 
authors, we avoided the inclusion of unnecessary 
climatic parameters that could place constraints on 
identifying climatically suitable habitat and chose 
a sub-set of seven key parameters based on expert 
opinion. These parameters were selected to represent: 
a) the most biologically meaningful predictors of 
koala distributions in terms of their thermoregulatory 
and nutritional needs (e.g. Moore and Foley 2000; 
Krockenberger 2003), and b) the most climatically 
influential requirements for Eucalyptus distribution, 
their primary food source (e.g. Box 1981; Hughes 
et al. 1996). These parameters were: Annual Mean 
Temperature; Maximum Temperature of the Warmest 
Period; Minimum Temperature of the Coldest Period; 


Annual Precipitation; Precipitation Seasonality 
(C of V); Precipitation of the Warmest Quarter; 
Precipitation of the Coldest Quarter (Table 1). 

Predictions: LGM range of P. cinereus 

It is indicated by Petherick et al. (2008) that the Last 
Glacial Maximum (LGM) in eastern Australia had two 
distinct cold and dry events (~ 30.8 k. cal. yr BP and 
21.7 k cal. yr BP) and was particularly severe. For 
example, temperatures were up to 8 ° C lower than 
today in mainland southeast Australia with summer 
precipitation reduced by at least 30% (Williams et 
al. 2009). In northeast Australia temperatures were 
3-7 ° C lower, with precipitation reduced to less 
than 50% of present values (Williams et al. 2009). 
We adjusted temperature and rainfall to generate 
two LGM predictions based on: 1) temperature 6 °C 
lower and rainfall 20% lower and 2) temperature 6 ° 
C lower and rainfall 40% lower. We chose regionally 
generalised climatic adjustments from approximate sea 
surface temperature (SST) reconstructions at the LGM 
(Turney et al. 2006; Van der Kaars et al. 2006; Barrows 
et al. 2007) and terrestrial and near-shore proxy data 
(Williams et al. 2009). We confined our predictions 
to the LGM due to the limitations in knowledge of 
the climate and the physiology of koalas at longer 
past time frames. We used the sub-set of seven key 
parameters to construct our LGM predictions. 

Results 

Predicted current range of P. cinereus 

The predictions under current climate conditions (Fig. 
3-a) identified a core range (10 th - 90 th percentiles) that 
falls within the modem distribution of koalas in the 
study area (Fig. 1). Some coastal areas were outside this 
core range, despite high densities of koalas currently 
being present in coastal areas such as southeast 
Queensland and eastern New South Wales. Conversely, 
the predicted entire range (1 - 100 th percentiles) included 
regions where there are currently no koalas such as 
Tasmania, Western and Central Australia. 

Predicted LGM range of P. cinereus 

With temperature reduced by 6 ° C and rainfall 
reduced by 20 % (Fig. 3-b), a highly restricted core 
range (10 th - 90 th percentiles) for koalas to southeast 
Queensland and coastal northern New South Wales 


Table I. Abbreviated Bioclimatic profile (percentiles) of Ph. cinereus for Queensland, New South Wales and Victoria 
showing the seven key parameters used for the Biomap LGM predictions. ('C ofV') = Coefficient of variation. 


Parameter 

Mean 

S.D 

5 

10 

25 

50 

75 

90 

95 

Max 

Min 

1. Annual Mean Temperature 

16.4 

3.39 

1 1.3 

12 

13.5 

16.5 

19.2 

21.2 

22.1 

24.2 

6 

2. Max.Temp. Warmest Period 

29 

3.67 

23.2 

24.1 

26 

28.8 

31.8 

34.4 

35.3 

37.7 

18.9 

3. Min.Temp. Coldest Period 

3.7 

2.65 

-0.6 

0.4 

2.1 

3.5 

5.3 

7.3 

CO 

oo 

15.7 

-4 

4. Annual Precipitation 

863 

327.58 

447 

493 

601 

796 

1061 

1333 

1528 

2480 

234 

5. Precipitation Seasonality (C ofV) 

36 

18.43 

10 

14 

23 

35 

47 

59 

70 

122 

6 

6. Precipitation Warmest Quarter 

290 

153.28 

88 

1 15 

178 

255 

376 

510 

599 

101 1 

61 

7. Precipitation Coldest Quarter 

172 

83.02 

49 

78 

107 

159 

227 

293 

329 

61 1 

17 
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C. 


Figure 3. Bioclim models of (a) modern predicted entire 
range and core range, (b) LGM predicted entire range 
and core range with temperature reduced by 6 ° C and 
rainfall reduced by 20%, (c) LGM predicted entire range 
and core range with temperature reduced by 6 ° C and 
rainfall reduced by 40%. 


was indicated. The entire predicted range (l-100 th 
percentiles) incorporated most of Queensland, some 
of coastal and northern New South Wales, small areas 
of coastal Victoria, South Australia and Western 
Australia and an area in the central Northern Territory. 
When the rainfall was further reduced by 40% (Fig. 
3-c), the predicted core range contracted away from 
the coast to highland areas in southeast Queensland 
and northern New South Wales. The entire predicted 
range contracted to eastern Queensland and New 
South Wales and highland areas of Victoria. 

Discussion 

Change in range boundaries 

We were able to model the potential distribution of P. 
cinereus at the LGM to demonstrate range contractions 
to areas that may have provided climatic refugia 
suitable for koalas during the severe conditions of past 
glacial periods. Our predictions concord with Korsten 
et al. (2009) whose mitochondrial DNA sequences 
for the brown bear Ursus ardos suggest that a small 
number of bears were restricted to a single refuge 
area during the last glacial maximum. The refugial 
hypothesis that populations retracted during glacial 
maxima is further discussed for wolves Cams lupus 
(Weckworth et al. 2010) and European deer Cervus 
elaphus and Capreolus capreolus (Sommer and Zachos 
2009). Phascolarctos cinereus fossils dated prior to the 
LGM are recorded in Western Australia’s Mammoth 
Cave (Glauret 1910) and Devil’s Lair (-31-43 thousand 
years ago) (Balme et al. 1978) as well as Madura 
Cave (>24-< 101 thousand years ago) (Lundelius and 
Turnbull 1982). We predicted a potential climatically- 
induced range contraction and eventual disappearance 
when we reduced rainfall by 40% in southwest Western 
Australia (Fig. 3c) that may indicate that as climatic 
conditions became more severe, there were no climatic 
refugia for koala populations to be sustained during 
the LGM. However any assumptions should be treated 
with caution due to uncertainty in fossil dating and 
the possibility of future palaeoecological information 
that may come to light. The question of whether 
Aboriginal hunting played a role in the disappearance 
of koalas from this region remains unresolved. 

Future distribution of koalas: climatic 
refugia 

The Pleistocene (2.6 million to 10,000 years 
ago), which incorporates the LGM, is considered 
an important period for assessing the impact of 
climate change on biodiversity due to its multiple 
oscillations between cold, dry glacial and warm, wet 
interglacial conditions. During this period, climatic 
fluctuations had profound effects on the distribution 
of vegetation (Johnson 2006), particularly impacting 
tree species and their distribution. Drought-induced 
deaths of eucalypts in Central Queensland are 
suggested to be more severe in species with roots 
confined to upper soil layers (Fensham and Fairfax 
2007) and forest dieback is predicted to increase 
worldwide with climate change (Cunningham et 
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al. 2010). It can therefore be expected that future 
koala distributions will be affected not only by 
the direct effects of a warming climate in terms of 
their thermoregulatory limitations but also by the 
indirect effects of food tree availability (Hughes 
et al. 1996; Allen and Breshears 1998) and foliar 
chemistry (Moore et al. 2005; DeGabriel et al. 2010), 
particularly in regions that become hotter and drier. 

To counter the impacts on koalas of rapid 
climate change, it may become necessary to 
consider conserving adequate and climatically 
suitable contemporary habitat, plus undertake 
novel conservation strategies such as assisted 
colonisation/migration (McLachlan et al. 2007; 
Hoegh-Guldberg et al. 2008). The latter option 
could involve moving koalas outside their current 
distribution, to climatically and ecologically 
suitable future climatic refugia. For example, the 
current range modelled here suggests koalas could 
survive in Tasmania and the alpine regions of 
Victoria, areas which are presently too cold and/or 
wet for their climate envelope. It is of interest that 
our predictions included Kangaroo Island (Fig. 3 
a, b), which currently supports a large introduced 
koala population (Fig. 1). Similarly, a small 


number of koalas translocated from Kangaroo 
Island currently inhabit Yanchep National Park 
in Western Australia. They are reportedly healthy 
and consuming Eucalyptus rudis, a species endemic 
to Western Australia (H. Beswick pers. comm.). 
However the problems that have developed with 
the introduction of koalas to Kangaroo Island 
(Duka and Masters 2005) demonstrate that assisted 
colonisation would firstly require precautionary 
prior field-based scientific experiments to test 
for the ecological and physiological viability and 
impacts of such measures. 

Our koala case study provides an historical context 
for the future conservation planning of this species in 
the wild and demonstrates the sensitivity of koalas 
to changes in climate, highlighting the importance 
of climatic refugia during the LGM. We suggest 
that historically, koalas have utilised such refugia 
during the harsher conditions of glacial maxima. This 
historical context and predictions of potential past 
range contractions to climatic refugia during harsh 
climatic conditions demonstrates that in planning 
for the future conservation of species threatened by 
rapid climate change, the identification of future 
climate change refugia will be a vital component. 
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